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Floquet theory is a powerful tool in the analysis of many physical phenomena. However, in 
its original formulation it is limited to linear systems with time-periodic coefficients. Here, we 
extend the theory by proving a theorem for the general class of systems including linear operators 
commuting with the period-shift operator. This class of systems ultimately includes those whose 
state depends on their full dynamical history. The present theorem greatly expands the range of 
applicability of Floquet theory to a multitude of phenomena that were previously inaccessible with 
this type of analysis. As an example, we provide the Floquet analysis of a Brownian particle with 
memory. 

PACS numbers: 05.45.-a, 02.30.Oz, 02.30.Sa, 84.30.Bv 



Floquet theory is a century-old theory that allows an 
easy analysis of the solutions of linear differential equa- 
tions that have periodic coefficients [TJ [5] . These types of 
equations are very common in several areas of science and 
technology, and as such the applications of Floquet the- 
ory range from quantum [3J H] to classical |S] physics, 
chemistry [7] , electronics [SJ [S] , dynamical systems [TU] 
and many more. However, its original formulation is lim- 
ited to linear systems with periodic coefficients. Conse- 
quently, it is a powerful tool to study nonlinear perturba- 
tions, noise, and stability of systems depending instanta- 
neously on time, and admitting periodic steady states. 

This leaves out many important phenomena that can- 
not be simply described using these types of equations. 
This is the case, for example, of all systems whose state 
depends on their past dynamics. These systems are the 
norm rather than the exception in natural phenomena 
[TT] and their applications span a wide range of problems 
(see, e.g., [T^] for recent examples). It would be thus of 
great value to extend Floquet theory to such systems. 

There have been some previous attempts to generalize 
Floquet theory to systems with memory. The most im- 
portant one is certainly the work of Hale and Verduyn 
Lunel [13] . They have derived a generalization of Floquet 
theory for linear equations of the type 

dz 

^=C v {z,t} (1) 

where C-p is a linear operator with the most important 
restriction that it is T-periodic, i.e., £-p{z,t + T} = 
£-p{z,t} for any t and any vector function z. However, 
this set of linear operators clearly does not cover the en- 
tire class of systems with memory. In fact, a general 
linear system with memory can be characterized by the 
linear operator 

C T {z,t}= f K(t;r)z(r) dr (2) 

Jt—r 



where r > and K is a matrix such that the integral 
of K(t; t) with respect to r is bounded for every t, and 
K(t + T; t + T) = K(f ; r). The latter condition (some- 
times called bi-periodicity) is not enough to guarantee 
the periodicity of Cj- for any z, but it is just true when z 
is T-periodic. However, by defining the time translation 
operator T as the linear operator T {v(t)} = v(t + T), it 
is simple to prove that Cj- commutes with T . This simple 
remark suggests that Floquet theory could be extended 
to systems including all linear operators commuting with 
T, thus allowing for the analysis of systems with memory. 

In this letter, we indeed state and demonstrate a the- 
orem that extends Floquet theory to such cases. The 
proof of the theorem relies only on well-known concepts 
of functional analysis, such as topological vector space 
and Banach algebra [13] • We then apply the theorem 
to the important problem of a Brownian particle with 
memory, where the latter is due to the presence of a fric- 
tion retardation function. In particular, we show how 
the Floquet theory enormously simplifies the study of 
the dynamical properties of this system. 

Let us then start by considering nonlinear dynamical 
systems described by a differential equation of the type 

^ = f(y,i)+£V{g(y,t),i} (3) 

where y : R — > R™, f , g : R™ +1 — > R" are nonlinear vector 
fields with f T-pcriodic with respect to t, and Of is 
a generic, linear integro-differential operator such that 
£V{g(y>*))*} is T-periodic if y(t) is T-periodic. We 
assume that ^ admits of a T-periodic solution ys (limit 
cycle). Dynamical systems with memory are included in 
([3]) as convolution operators 

The study of the perturbation of the limit cycle ys is 
important both for assessing the stability of the solution 
and for analyzing the effect of small-amplitude forcing 
terms such as Langevin sources used to represent fluctu- 
ation effects. Perturbations are studied by linearizing pi) 
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around ys yielding a linear periodic time- varying system 



dz 

di. 



A(i)z + £ r {z,t} + b(t) 



(4) 



where z : R — > R n , A(t) is the T-periodic Jacobian ma- 
trix of f(-,t) evaluated in the limit cycle, and h(t) is the 
perturbation forcing term. Finally, Cj- is the linear oper- 
ator defined as Cj- {•, t} = Cj- {B(i)-}, where B(i) is the 
Jacobian matrix of g(-,t) evaluated in ys- The solution 
of the homogeneous equation associated to Q (i.e., the 
case b(f) = 0) is denoted as zo(t). The only limitation we 
assume for the linear operator in Q is that it commutes 
with the time translation operator T, and it acts on func- 
tions z defined at least on the interval [— r+s, s+T] where 
s is the initial time instant for Q. 

Denoting as £(T) the set of a// the linear operators Cj- 
that commute with T, a trivial verification of the defini- 
tion conditions shows that £(T") is a Banach algebra gen- 
erated by T on the topological vector space H, a Hilbert 
space made of functions v(t) : [— r + s, s + T] — s- R n . 

After this brief but necessary introduction we can then 
formulate the main result of this contribution in the form 
of the following theorem: 

Theorem (Generalized Floquet theorem). Let us con- 
sider the homogeneous part of system Q with Cj- £ 
£(7~), whose space of solutions is spanned by the state 
transition matrix X(t; s) (for s < t < s + T). The state 
transition matrix has then size n x p with p < +oo, 
and can be written as X(t;s) = M(t; s) cxp[F(i — a)], 
where M(t; s) <E R nxp is T-periodic in both variables (i.e., 
M(t + T; s) = M(i; s) and M(t; s + T) = M(t; s) J and 
F G |RP x p is « constant matrix. 

Proof. The first part of the proof is devoted to the 
definition of the size of X(t; s). Let us consider the op- 
erator 



{A(t)-}-£r {.,*}. 



(5) 



£'V is a linear operator acting on elements z £ H, the 
Hilbert space spanned by the functions z : [—r + s, s + 
T] —> R n . Since the functions z are defined over the 
finite interval [— r + s, s + T], % has an infinite countable 
dimension and according to [2], ker {£"7-} is a countable 
Hilbert space of dimension p < +00. As consequence 
X(t;s) € R nxp . 

Let us now turn to the form of X(t; s). We immediately 
see that £'V{X(t; s )} = because it is a state transition 
matrix for @, and clearly T{£"r{X(£; s)}} = T{0} = 
0. Since T commutes with £7-, with d Q /d£ a for any 
a £ N, and, as can be easily proved, also commutes with 
all the time periodic matrices we get 

T{£" r {X(t; s)}} = £" r {r{X(<; s)}} = 0. (6) 

Equation ^ implies that T{X(t;s)} is a state tran- 
sition matrix of Q. Being formed by the basis elements 



of the kernel of the linear operator £"tj the state tran- 
sition matrix 7~{X(i; s)} can be expressed as a combina- 
tion of the elements kj (t) of a base of ker{£"7-}, i-e. as a 
combination of the columns of X(£; s). From elementary 
representation theory, a representation of the linear oper- 
ator T is a matrix, therefore such combination is linear. 
These remarks lead to 



T {X(t; a)} = X(t + T: s) = X(t; a)C. 



(7) 



where the matrix C € IR pxp is constant because T com- 
mutes with d a /dt a . 

Since the exponential function is never null, without 
any loss of generality we can require that the form of 
X(i; s) is X(t; s) = M(t; s) exp[F(i - s)], where F e Rp x p 
is a constant matrix. Using this in ([7]) 

T {M(f ; s) e F ^ - s ) } = M(t + T;s)e F ( t + T - s ) 

= M(t; s) e F (* " s ) C (8) 

that is satisfied only if C = exp(FT) and M(t + T; s) = 
M(t; s), that is the result we set to prove. Repeating the 
proof with a time translation operator acting on s, we 
finally find that M(t; s) is periodic also with respect to 
s. This completes the proof of the theorem. 

Discussion. Because of the generalized Floquet theo- 
rem, a solution of the homogeneous part of Q is 



z (t) = M(t; 0) e F< = M(t; 0)Q e ul c (9) 

where D £ Rp x p is a diagonal matrix corresponding to 
the solution of the eigenvalue problem FQ = QD, and 
Cq = Qcq £ R ,n is a constant vector. Choosing Co = ej 
(the j-th element of the canonical unit vector base of R p ), 
Eq. ([9| yields the general form of the solution 



z (t)=T j (t)e X i t 



J 



,p, 



(10) 



where Xj is the j-th diagonal element of D and Yj(t) 
is a T-periodic vector function. Dropping the index j, 
substituting ( 10 ) into Q yields 

^ + Ar = A(t)r + c- Ai £ r {e A< r,t}, (11) 



where the last term on the right hand side is T-periodic 
because all the other terms are. The asymptotic stability 
analysis of the limit cycle ys, therefore, clearly depends 
on the computation of all the eigenvalues of F: this pro- 
cedure depends on the explicit form of £7-, and in general 
leads to a set of p Floquet exponents A (notice that ac- 
cording to the generalized Floquet theorem p may also 
be infinite). 

Another important remark is that direct substitution 
proves that the general solution of Q can be written as 



«(t)=X(t;0)co+ / X(t;r)b'(r) dr 
Jo 



(12) 
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where Co is a constant vector of size m, and h'(t) is a 
solution of 

X(0; O)b'(t) = b(t) + C r I J X(t; r)b'(r) dr J 

£ T {X(t;r)}b'(r)dr. (13) 

Jo 

Finally, the results of the theorem can be extended to 
operators of type ^ for r — > +00 under the following 
considerations: when r — > +00, because of the properties 
of K, K(<, r) — > for r — > — 00 and any t. For any e > 

there is an r = r(e) such that /_ r K(t; r)yg dr < e 

being ys T-periodic and integrable. Thus, for any e we 
can approximate a infinite memory system with ([2|. In 
this sense our theorem guarantees that the solution of 
the homogeneous part of system Q is given by ( [To| for 
t G [— r + s, s + T] and for any r < +00. 

Example. As an important example of application we 
consider the analysis of a driven Brownian particle with 
memory, where the memory part is due to the presence 
of a friction retardation function as discussed in |15j . 
The particle is subject to an external field with poten- 
tial energy U(x). According to [15], the memory particle 
Langevin equation is 



dv 

"dT 



f T(t, r)v(r) dr - V x f/(x) + (14) 

Jt-r 



where v = dx/di, is the zero average Langevin 

source, and for a particle in thermal equilibrium 
the friction retardation function is |15j T(t, r) = 
^^(£ t (£)£(t)), where (•) is the ensemble average, T 
the transpose operation, m the particle mass, 6 the ab- 
solute temperature, &b the Boltzmann constant, and 
r > such that f_ r T(t, t)v(t) dr ps 0. Consider- 
ing an Ornstein-Uhlenbeck process where (£ T (i)^(7-)) = 
kkBd fi exp(— k\t — r|) [16] and the mobility is taken veloc- 
ity dependent, i.e. \i = ^[(v}(t)], the friction retardation 
function reads 



m 



-k\t 



(15) 



where k = l/r n , r n is the noise correlation time, and v 



the magnitude of v. Taking the average of ( 14 1 and using 



(15 



d(v 



k 

dt m Jt _ r 

(V x tf(x)) 



M [(«)(r)]e- fc l f - r l (v)(r) dr 



(16) 



Equation ( |16| is coupled to (v) = d(x) /dt thus leading to 
a problem of type ([3]) with unknown y = [(v) T (x) T ] T . 

Linearizing around the periodic solution for £ = 0, the 
homogeneous variational memory equation reads (since 



T <t) 

d6v 

~dt~ 



k 
m 



/ ii'(T)e- k{ t- T ) 5v{t) dr-G(«)5x 

Jt-r 

(17) 

where ,/(i) = m[(v s )(<)] + V vA i[(v s ) (f)] T (v s )(t) and 
G(*) = V x £/((x s )). Clearly, ([l7|> is completed by 5\ = 
d<5x/d£ therefore z = [5v T <5x r ] T . According to the 
generalized Floquet theorem and the previous discussion, 
the solution of ( [l7] ) is in the form (10). Decomposing 
r(t) = [rj rJ] T into the velocity and position compo- 
nents, and being p! and T-periodic, the explicit form 



of ( 11 ) for the v part reads 



dv v 
dt 



Ar„ 



-G(i)r a 



£ 



kfl h r V j-h 1 - e 



to(/c + A + icjj) 



c i£ -V (18) 



where i = \/— 1, j G T, ujj = ju>o (ujq — 2%/T), and 
with c\j we indicate the j'-th harmonic amplitude (in ex- 
ponential form) of the generic T-periodic function a(t). 
Eq. |l8| is completed by dr x /dt + Xr x = r„. Tak- 
ing r large enough, for all A satisfying k > — Re(A), 
exp[— (k + A + iu)j)r] can be neglected. It is worth notic- 
ing that the last assumption is always satisfied for all 
A with real part close to or larger, and consequently 
responsible for bifurcations and instabilities. So, by 
looking at this subset of critical exponents, we neglect 
exp[— (fc + A + iwj)r], expand G(t) in Fourier series since 
it is also T-periodic, and finally by balancing the har- 



monic components in ( 18 ) we find the nonlinear eigen- 
value problem 



(iwj + A) r v j + 



m(k + A + vjjj) 



t l h r v,j-h 



^ Ghr Xi j-h — 0. 



(19) 



Of course, ( 19 1 is completed by the frequency domain 



version of the linear relation linking r x (t) and r v (t): 



Multiplying the infinite set ( 19 ) by 



k + A + icjj we obtain a polynomial eigenvalue problem 
(PEP) of order 2. The r order PEP has in general r x n 
(with n the number of equations) eigenvalues [17l [18], 
meaning that in our case p = 2(n/2) + n/2 = 3n/2. In 
order to numerically tackle this problem, we first reduce 
the PEP to a standard eigenvalue problem with the meth- 
ods proposed in [TTJ [TS] and then calculate the harmonics 
fj as in [T9] . 

Specifically, we consider the two-dimensional Brown- 
ian particle system studied in [50], and modify it by 
including memory. The potential energy is f/(x) = 
0.5o) 2 x where x\ and X2 are the 



§.h(<jj\x\ + U)2 X 2 



2_2-\ 
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FIG. 1. Arnold tongue in the (o>i/o>2, a) plane for the driven 
Brownian particle with memory. Parameters: /3 = 1, oh. = 2. 



two components of the position x of the particle, ui\ and 
UJ2 are two constant parameters, and u) = diag {cJi, 0T2}. 
This means that G(t) = u) 2 is constant, thus reducing 
the second series in |T9| ) to ui 2 r x j. According to [20], 
we use a corrected Rayleigh type model for the particle 
viscosity: 



H[v{t,T)]=j R [v(T)}5(t-T) 1 -+g' 



(20) 



where Jr(v) — —(a — f3v 2 ), a, (3 > and g' are param- 
eters. We also define the parameter g — g'/k; notice 
that for constant <?' , a memoryless system is obtained for 
k -> + 



This means that ( 19 ) takes the form 



+ A + 



k + A + iujj 



ft 

.2- 







(21) 



where 7^(t) = 7R.[v s (i)] + V v 7R[vsW] T vs(i). Since 
n = 4, according to the previous discussion p = 6. We 
implement a numerical solution of the second-order poly- 
nomial eigenvalue problem using N = 30, thus expanding 
(21| into a system for j = —N,...,N. Furthermore, 



the limit cycle solution is calculated in the frequency 
domain exploiting the Harmonic Balance [T9] numerical 
technique, again making use of 30 harmonics. 

Concerning the numerical results, with respect to the 
memoryless case a stable equilibrium (originally unsta- 
ble) is found in the origin of the phase space. A bi- 
furcation diagram in the parameter space (0)1/0)2,0:) is 
plotted in Fig. [T] for k = 7 and g = 4. An Arnold tongue 
is shown, separating a stable equilibrium for the lower 
a values (not present in |20j ) from two symmetric, sta- 
ble limit cycles (the inner part of the Arnold tongue). 
Out of the Arnold tongue, two similar strange attractors 
appear generated by the collapse of the two limit cycles 
as discussed in [2U]. The bifurcation curves correspond 
to fold bifurcations, i.e., a Floquet exponent with zero 
imaginary part crossing the zero real part boundary. 




(£>J(£> 2 

FIG. 2. Arnold tongue in the (0)1/0)2, ct) plane for the driven 
Brownian particle with memory as a function of g; parame- 
ters: P = 1, 0)1 = 2, k = 7. Inset: Bifurcation curve in the 
triple point (0)1/0)2 = 1) in the (fc, a) plane as a function of 
g'; parameters: /3 = 1, 0)1 = 2. 



The evolution of the Arnold tongue as a function of g 
and k, and the bifurcation curves associated to the triple 
point 0)1/0)2 = 1 as a function of g' , are shown in Fig. [2] 
The curves tend to zero for k — > +00 since in this limit 
the memory part of the system disappears. 

Conclusions. In conclusion, we have stated and 
proved a theorem that generalizes the standard Floquet 
theory to systems whose state depends on the full past 
dynamics. As an example of application we have ana- 
lyzed the dynamical properties of a Brownian particle 
with memory. Such a theorem greatly expands Floquet 
analysis to systems that were previously inaccessible by 
such theory, making it a tool with even greater poten- 
tial in several areas of science and technology. Also, it 
is worth mentioning that the theorem we have proved 
applies also to Hamiltonian systems with non-local (al- 
beit spatially periodic) interactions. Indeed, the stan- 
dard Bloch theorem of solid-state physics [5T] is a simple 
application of Floquet theory, obtained by replacing the 
time variable with the space variables. By the same to- 
ken, our theorem is a generalization of Bloch theorem 
to non-local spatial correlations which can be generally 
cast in the form of Eq. ^ , by replacing the time variable 
with the space variables. We thus expect it to be of great 
value to study the equilibrium properties of systems that 
are quite common in condensed matter physics. 
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